Energetics of superconductivity in the two dimensional Hubbard model 
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The energetics of the interplay between superconductivity and the pseudogap in high temperature 
superconductivity is examined using the eight-site dynamical cluster approximation to the two 
dimensional Hubbard model. Two regimes of superconductivity are found: a weak coupling/large 
doping regime in which the onset of superconductivity causes a reduction in potential energy and 
an increase in kinetic energy, and a strong coupling regime in which superconductivity is associated 
with an increase in potential energy and decrease in kinetic energy. The crossover between the two 
regimes is found to coincide with the boundary of the normal state pseudogap, providing further 
evidence of the unconventional nature of superconductivity in the pseudogap regime. However the 
absence, in the strongly correlated but non-superconducting state, of discernibly nonlinear response 
to an applied pairing field, suggests that resonating valence bond physics is not the origin of the 
kinetic-energy driven superconductivity. 
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The high transition temperature superconductivity ex- 
hibited by layered copper-oxide materials has been an 
important topic in condensed matter physics since its 
discovery in 1986^ Broadly speaking, two views are cur- 
rently held. One is that despite the various anomalous 
features of the materials the superconductivity may be 
understood in more or less conventional Bardccn-Coopcr- 
Schrieffer (BCS) terms as arising from the exchange of a 
pairing ('glue') particle, most likely of magnetic origin3 
An alternative view is that the superconductivity is an 
intrinsic property of a strongly correlated state of mat- 
ter that should not be interpreted as arising from the 
exchange of a well-defined excitation^ 

The issue may be cast in energetic terms. In the con- 
ventional BCS view, the driving force for superconduc- 
tivity is in essence a reduction of potential energy: by 
forming the superconducting state the electrons can take 
greater advantage of an attractive term in an interparticle 
potential. Changing the wave function to reduce the po- 
tential energy however costs kinetic energy, so that in the 
weak coupling limit the change from normal to supercon- 
ducting states leads to an increase in the kinetic energy.— 
In an alternative view^ the driving force for supercon- 
ductivity is an optimization of kinetic energy: by forming 
the superconducting state the electrons can move more 
easily through the crystal despite their need to avoid the 
other electrons. In this case going from the normal to 
the superconducting state lowers the kinetic energy and 
one expects that the potential energy increases. 

The repulsive-U Hubbard model on the two dimen- 
sional square lattice is widely believed£££ to contain the 
essential physics of high-T c copper-oxide superconductiv- 
ity. It is defined by the Hamiltonian 
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FIG. 1. (Color online) Phase diagram of two dimensional 
square lattice Hubbard model in plane of density n and 
interaction strength U/t at inverse temperature (3 = 60/i 
as obtained in 8-site cluster dynamical mean field the- 
ory. Mott insulator at half filling for U/t ^ 6.4 indi- 
cated by heavy bar (green online); superconducting region 
indicated by circles (black online), pseudogapped but non- 
superconducting region, diamonds (blue online) and Fermi 
liquid non-superconducting state by squares (yellow online). 
Boundary of normal state pseudogap, defined as in Ref. HQ . 
indicated as dashed line (purple online). Trajectories along 
which the energy is computed are shown as arrows. 



repulsive (U > 0) version of the model has been shown 
rigorously to have a d x 2_ y 2 superconducting ground state 
in at least some regions of the U, n phase diagram.— ~— 

In this paper we investigate the electronic energy E = 
(H), decomposed into kinetic K and potential V terms 
as E = K + V with 

K=J2^k{c{ a c ka ) = 2T^(e fe - M )Tr [r 3 G(fc, u n )] (2) 



kcr 



Here i labels the sites in a lattice and k a momentum in 
the corresponding Brillouin zone. The two dimensional 



V=U^(n itnii ) = 2Ty^Tr[B(k,u; n )G(k,oj n )] (3) 
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In the second equality we have used standard formulae to 
reexpress the expectation values in terms of the Nambu 
matrix Matsubara frequency electron Green function G 
and self energy S. 

The energetics of superconductivity have been pre- 
viously studied. One important class of approaches 
has used variational wavefunctions, often starting from 
a single Slater determinant with doubly occupied sites 
then being projected out.— ~— Information about pair- 
ing comes from comparing results obtained from free 
fermion and BCS-paired starting points. These works in- 
dicated that pairing was present for dopings from x = 
to x ~ 0.25 and that over most of the phase diagram the 
kinetic energy of the paired state was lower than that of 
the unpaired state. However, variational results are con- 
strained by the choice of variational space; in particular 
by the choice of projective BCS-type wave functions. 

Another important class of theoretical approaches in- 
volves phenomenological spin-fermion models.— ~— In 
these approaches it is assumed that the important physics 
arises from the interaction of electrons with spin fluctua- 
tions (treated as bosons but with boson self-energy effects 
arising from coupling to fermions playing a crucial role). 
These models are amenable to semi-analytic treatment. 
Their analysis revealed that in the strong coupling limit 
the superconducting state could have lower kinetic en- 
ergy than the normal state. However, these models do 
not fully capture the strong correlation effects associated 
with the Mott transition or the formation of the pseudo- 
gap, and rely on assumptions about the most physically 
relevant interactions. 

We use the dynamical cluster approximation (DCA) 
version of dynamical mean field theory^ to evaluate 
Eqs. |2I3| for the two dimensional rcpulsive-U Hubbard 
model with Ek = —2t(cosk x + cosk y ). In the DCA the 
Brillouin zone is tiled with TV patches and the electron 
self energy is taken to be piecewise constant, with a differ- 
ent value in each sector of momentum space. The sector 
self energies are obtained from the solution of an auxil- 
iary quantum impurity model with parameters fixed by 
the Hubbard interaction and a self-consistency condition 
discussed in detail in Ref.[26l The method yields a d x 2_ y 2 
superconducting state^l - — For the Hubbard model the 
method becomes exact as N — > oo and considerable ev- 
idence is now available^ - — concerning the status of the 
finite TV results achievable numerically. Here we study 
the case TV = 8, which has been shown to be large enough 
for the results to be representative of the infinite cluster 
size Hmil£ but small enough to enable calculations of the 
necessary accuracy4£ 

We obtained the superconducting kinetic and potential 
energies KEg and PE$ from superconducting solutions 
obtained as described in Ref. and the normal state en- 
ergies KEn and PEn by solving the DMFT equations in 
the paramagnetic phase with the same code but subject 
to the constraint that the anomalous ((cc)) terms in the 
Green function and self energy vanished. Our results are 
obtained using the CT-AUX version 4 ^ of the continuous- 



time quantum Monte Carlo method^ 4 - with submatrix 
updates 4 ^ and an extension to superconductivity^ The 
energy differences are found to be very small and careful 
attention to the high frequency behavior is required for 
reliable results. The submatrix methods are essential in 
obtaining data of the requisite accuracy. 

Fig. Q] shows the phase diagram obtained from the 
TV = 8 DCA method in the interaction strength and 
doping planed along with two arrows indicating the 
parameter-space trajectories along which energies are 
computed in this paper. At U > 6.4i and carrier con- 
centration n = 1 per site the approximation yields a 
paramagnetic ('Mott') insulating state which is at lower 
temperatures unstable to antiferromagnetism. As elec- 
trons are removed the state evolves to a conventional 
Fermi liquid metal via an intermediate 'strange metal' 
phase characterized by a 'pseudogap', a suppression of 
electronic density of states in the (0, ir) region of the 
Brillouin zonei 10 ' 40 ' 46 " — Superconductivity is found in a 
stripy near to the Mott insulator but separated from 
it by a region of pseudogapped but nonsuperconducting 
states^ At carrier concentration n = 1 (vertical arrow) 
the ground state of the model is believed to be antifer- 
romagnetic at all U. The n = 1 results were obtained 
by suppressing long-ranged antifcrromagnetic order (al- 
though short-ranged antifcrromagnetic correlations are 
still present) and are representative of the properties of a 
metastable state. They are included because the qualita- 
tive properties are seen to be the same as in the doping- 
driven transition but the particle-hole symmetry at n = 1 
permits the acquisition of much higher quality data, en- 
abling a clearer view of the phenomena. 

The two panels of Fig. [5] show the energy differences 
obtained by subtracting the superconducting and normal 
state energies computed at inverse temperature /3 = 60 /t 
along the two parameter-space trajectories shown by the 
arrows in Fig. [TJ i.e. crossing the superconducting re- 
gion by varying the interaction strength or varying the 
carrier concentration. The results obtained along the two 
trajectories are remarkably similar, although the absence 
of a fermion sign problem at n = I means we arc able 
to obtain much better statistics in this case. The con- 
densation energy is of order O.OOlt, although the changes 
in kinetic and potential energy separately are typically 
much larger, especially in the pseudogap regime. 

For interactions or carrier concentrations smaller than 
required to produce a normal state pseudogap ; 10 ' 40 ' 49 the 
energetics are consistent with the standard expectations 
of weak coupling superconductivity: as the material en- 
ters the superconducting state the potential energy de- 
creases and the kinetic energy increases. The boundary 
of the normal state pseudogap marks a significant change 
in the energetics of superconductivity: once the pseudo- 
gap regime is entered, the kinetic energy decreases and 
the potential energy increases on entering the supercon- 
ducting state. Further, inside the pseudogap regime the 
superconducting/normal changes in potential and kinetic 
energy become much larger in magnitude, showing that 
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FIG. 2. (Color online) Differences in total, kinetic and poten- 
tial energies (per site, in units of hopping t) between normal 
and superconducting states, obtained as described in the text 
at density n = 1 varying interaction strength (upper panel) 
and as function of density at fixed interaction strength U — 6t 
(lower panel). 



the onset of superconductivity leads to a significant re- 
organization of the energetics of the pseudogap states. 
The change in character of the superconductivity at the 
pseudogap line is consistent with the finding of Yang 
et al— that the superconductivity exists in a dome with 
the maximal transition temperature occurring where the 
superconducting and pseudogap phase boundaries inter- 
sect. 

Our results differ from previous dynamical mean field 
analyses. Ref. HH (N = 4 study of the Hubbard model) 
and Ref. [H (N = 4 study of the t- J model, with an addi- 
tional 'EDMFT' approximation) found that most of the 
energy gain on entering the superconducting state came 
from changes in the interaction term, although Ref. [56| 
found that the behavior of the kinetic energy was differ- 
ent at large than at small doping. Three possible origins 
for the discrepancy are the use of the non-crossing ap- 
proximation ("NCA") impurity solver in Refs. l55l andl56l 
rather than the numerically exact CT-QMC method, the 
use of the N = 4 approximation, rather than the N = 8 
approximation used here, and the study of the t-J rather 
than Hubbard model in Ref. [56|. Singh^ has questioned 
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FIG. 3. (Color online) Anomalous expectation value in sector 
K = (0, 7r) plotted against pairing field r\ at doping x = for 
interaction strengths indicated. 



the relevance of computations based on the t-J model, be- 
cause of apparent violations of the virial theorem which 
may be traced back to the fact that some parts of the 
electron kinetic energy are included in the 'J' coupling. 

The 'potential energy-driven' nature of the supercon- 
ductivity found at larger dopings and at weak couplings 
is consistent with the notion that in these regimes the 
superconductivity is relatively conventional. The change 
in energetics as the pseudogap boundary is crossed sug- 
gests that at stronger couplings or lower dopings the 
superconductivity becomes unconventional. One in- 
fluential model of unconventional superconductivity is 
the resonating valence bond (RVB) idea of Anderson^ 
which was motivated in part by the possibility that the 
physics of the cupratcs could be understood in terms 
of a very strong coupling limit of the Hubbard model. 
There, configurations with two electrons on a site could 
be projected out so that the only important term in the 
energy was the kinetic energy term and superconductiv- 
ity (and indeed all other interesting physics) is necessarily 
'kinetic energy driven'. 

Anderson's original RVB idea, and subsequent recog- 
nition of an SU (2) symmetry in the RVB wave func- 
tion which might be weakly broken by doping or small 
finite U corrections^ implied that superconducting cor- 
relations were present (but not active) in the Mott insu- 
lating state and the strongly correlated but not supercon- 
ducting state which separates the superconducting and 
insulating regimes in Fig. [TJ To test this hypothesis we 
applied a pairing field t/k {c ! K ^c ! K ^ + ck^cki) in our calcu- 
lations and computed the effect on the superconducting 
order parameter (ck^ck\)- We expect that if a near- 
SU(2) symmetry existed, then applying a small pairing 
field to a state which is non-superconducting but is near 
the phase boundary would provide a rapid increase in 
the pairing amplitude, which would then saturate to a 
value characteristic of the superconducting state. Fig. [3] 
shows that this is not the case. On the weak coupling 
side (U — 4.2), applying a pairing field leads to the be- 
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FIG. 4. (Color online) Temperature dependence of kinetic 
energy at n = 1 for U = 5.0 (lower panel) and U = 5.8 
(upper panel) in normal (filled squares, red dashed line) and 
superconducting state (open circles, black solid line). 



havior expected near a second order phase transition: a 
rapid increase in (cc) reflecting the enhanced susceptibil- 
ity, followed by a saturation to values similar to those 
found in the nearby superconducting state. However, on 
the strong coupling side the situation is different. Just 
at the phase boundary U = 6.0 the situation is similar to 
that found at weak coupling, but for any larger U the (cc) 
vs r/ curve is linear with small, weakly U-dependent slope. 
The similarity of the U = 6.2 and 6.4 results, and the dif- 
ference of both of these to the U = 4.2 trace, indicates 
that precursor effects arc very weak as the superconduct- 
ing phase is approached from the pseudogap indicating 
that the pseudogap state has no strong tendency towards 
superconductivity. We infer from this calculation that 
the origin of the kinetic energy-driven behavior is not a 
signature of pairing correlations pre-existing in the wave 
function. 

It is interesting to consider the normal - superconduct- 
ing energy differences in the context of the energetics of 
the pseudogap state itself. The two panels of Fig. [4] show 
the temperature dependence of the kinetic energy com- 
puted for a relatively weak coupling, U = 5. Of, (lower 
panel) and relatively strong coupling, U = 5.8t, (upper 
panel) . We see that in the weak coupling case, the kinetic 
energy decreases as the temperature is lowered, and the 
onset of superconductivity reverses this decrease, while 
in the strongly coupled case the kinetic energy increases 
as temperature is lowered but the onset of superconduc- 
tivity again reverses the temperature dependence. 

We finally consider the observability of these effects. 
Norman and co-workers noted that the difference be- 
tween normal-state and superconducting state photoe- 
mission spectra could be analyzed to obtain estimates of 
the normal-superconducting change both in total and in 
kinetic energy^ although the analysis is complicated by 
the need to accurately monitor small changes occurring 
over wide energy and momentum ranges. Specific heat 
data can also be used to infer the condensation energy, 



although care must be taken both to extrapolate the nor- 
mal state to temperatures less than the transition tem- 
perature and to include fluctuation effects^ In narrow- 
band systems such as cuprates, an approximate relation 
between the kinetic energy and the frequency integral 
of the optical conductivity exists^ir— Interestingly the 
idea of examining superconductivity-induced changes in 
the optical integral seems to have entered the high-T c lit- 
erature first in the context of the 'hole superconductiv- 
ity' model of Hirsc h 64 ' 65 and then by Chakravarty and 
collaborators^ in relation to the interlayer coherence 
mechanism of Anderson^! These works motivated sev- 
eral experimental groups to examine changes in optical 
conductivity across the normal-superconducting phase 
boundary^r— Unfortunately, the value of the optical 
spectral weight depends on the frequency up to which 
the conductivity is integrated, and the appropriate up- 
per cutoff may be different in the normal and supercon- 
ducting state j 23 ' 24 also even if the conduction band con- 
tribution to the optical sum rule could be determined, 
the relation between this and the kinetic energy is only 
approximate, and the errors involved in the approxima- 
tion may be different in the normal and superconducting 
stated As can be seen from Fig. [4] the temperature de- 
pendent changes are only on the 1% level. 

In summary, our results indicate that the nature of 
the superconductivity depends crucially on location in 
the phase diagram. For the high doping/weak correla- 
tion side of the superconducting region, the energetics 
of superconductivity appear essentially conventional: on 
transition to the superconducting state the potential en- 
ergy decreases and the kinetic energy increases. On the 
low doping/strong correlation side of the superconduct- 
ing region, the energetics appears unconventional: on en- 
tering the superconducting state the kinetic energy de- 
creases and the potential energy increases. However, we 
do not find any indication that the nonsuperconducting 
pseudogap state has any significant pairing correlations, 
casting doubt on an RVB interpretation of the pseudo- 
gap. Interestingly, the crossover between the two regimes 
occurs essentially at the point at which the pseudogap 
becomes visible in normal state quantities and also in- 
terestingly the changes in the individual components of 
the energy become much larger in this unconventional 
regime, suggesting that superconductivity causes a sub- 
stantial rearrangement of the pseudogap electronic state. 
In other words, the superconductivity and pseudogap are 
competing phases. 
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